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Abstract — Due to the very high nonlinearity of the power system, traditional analytical methods take a lot 
of time to solve, causing delay in decision-making. Therefore, quickly detecting power system instability 
helps the control system to make timely decisions become the key factor to ensure stable operation of the 
power system. Power system stability identification encounters large data set size problem. The need is to 
select representative variables as input variables for the identifier. This paper proposes to apply wrapper 
method to select variables. In which, Binary Particle Swarm Optimization (BPSO) algorithm combines 
with K-NN (K=l) identifier to search for good set of variables. It is named BPSO&l-NN. Test results on 
IEEE 39-bus diagram show that the proposed method achieves the goal of reducing variables with high 
accuracy. 
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I. INTRODUCTION 

Modem power systems suffer from operating pressure 
very close to a stable boundary limit, while power systems 
are always faced with abnormal agitation that can easily 
cause instability (Makarov, Reshetov, Stroev, & Voropai, 
2005). Due to the very high nonlinearity of the power 
system, traditional analytical methods take a lot of time to 
solve, causing delay in decision-making. Therefore, 
quickly detecting power system instability helps the 
control system to make timely decisions become the key 
factor to ensure stable operation of the power system. 

In recent years, the identification method has been applied 
as an alternative to solving difficult problems that 
traditional methods of analysis cannot solve in terms of 
computational (Dong, Rui, & Xu, 2012; Mitra, Benidris, 
Nguyen, & Deb, 2017; Oliveira, Vieira, Bezerra, Martins, 
& Rodrigues, 2017; Zhang, Xu, Dong, Meng, & Xu, 
2011). By learning the database, the nonlinear input/output 
relationship between the electrical system operating 
parameters and stability can be quickly calculated. 
However, if the identifier works fast, the number of input 
variables must be characteristic variables. Therefore, the 
selected variables must be characteristic, eliminate 
redundant variables and cause noise. 

The selection of variables of previously published works 
mainly applies the filter method that has to go through 
many stages (Dong et al., 2012; Kalyani & Swarup, 2013; 
Swarup, 2008). With this method, the first is based on the 
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statistical criteria to filter out the good variables first, then 
ask the identifier to evaluate the accuracy to find the 
appropriate set of variables. This article recommends a 
wrapper method to select variables. In which, BPSO 
(Binary Particle Swarm Optimization) algorithm combined 
with K-NN (K=l) identifier to perform the process of 
searching for good set of variables. It is named BPSO&l- 
NN. This helps to reduce the step of finding the variable 
set compared to the filter method and improves the 
automation of the process of finding the variable set. Test 
results on IEEE 39-bus diagram show that the proposed 
method achieves the goal of variable reduction with high 
accuracy. 

H. VARIABLE SELECTION 

In the design stages of a stable diagnostic model of the 
power system using the identification method, the 
selection of characteristic variables is very important 
because it directly affects the diagnostic accuracy of the 
model. This also helps to reduce the number of 
measurement sensors, reducing computation time for the 
model. This step defines a specific set of variables that 
represent the database for training. These initial 
characteristic variables are the input variable representing 
the operating parameters of the power system and covering 
the operating status of the power system. The 
characteristic variable of the power system in transient 
mode or dynamic mode is the change in generator 
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capacity, change of load capacity, change of power on 
transmission lines, voltage drop at nodes, ... right at the 
time of the incident. The output variable is labeled binary. 
‘0’ is unstable and ‘1’ is stable. 

Filter method: 

The Filter method applies the search algorithm with the 
statistical standard function such as Fisher, Divergence, ... 
to guide the search strategy for the best set of variables to 
achieve the target value (J) (Dong et al., 2012; Kalyani & 
Swarup, 2013; Swarup, 2008). The process of searching 
for the set of variables by the filter method is presented as 
Figure 1. This method shows that the result is a set of 
variables meeting statistical standards but does not show 
the identification accuracy and the number of necessary 
variables. Therefore, the filter method needs to add an 
identification accuracy assessment to select the last set of 
variables. 
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-► 
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Fig. 1: Diagram of Filter method 
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Fig. 2: Diagram of Wrapper method 


The Wrapper method is presented as in Figure 2. However, 
this method differs from the Filter method in that the target 
function considering the identity accuracy. Thus, this 
method results in a set of variables to be selected with a 
specific identification accuracy. 


The fitness function for selecting variables is presented in 
formulas (1), (2), and (3). 


ErrsF 

Fitness =- 

TotalF - sF 

(1) 

Acc S F(%) = A « S “'”P“ (%) 
TotalSample 

(2) 

ErrsF (%) = 1 — AccsF 

(3) 

where: 


Finess is the objective function. 



AccsF is the identification accuracy of the selected set 
of variables. 

AccSample is the number of the correct identifier. 
TotalSample is the total number of samples. 

TotalF: the total number of variables in the data set. 
sF is the number of selected variables. 

ErrsF is the identification error of the selected variable. 


III. BPSO ALGORITHM IN VARIABLE 
SELECTION 


PSO (Particle Swarm Optimization) is the optimal search 
algorithm proposed by Kennedy and Eberhart (Kennedy & 
Eberhart, 1995). In the PSO, each candidate answer of the 
problem is encoded as an instance that moves through the 
search space. The whole herd seeks the optimal solution by 
updating the position of each individual based on their 
own experience and on neighboring individuals. 

Generally, the vector Xi = (xu, xn, ..., x,d) used in the PSO 
represents the position of the i th instance. The vector v, = 
(vu, Vi 2 ,..., Vm) is used in the PSO to represent the velocity 
of the i th instance. D is the size of the search space. 
During the search, the best position for each previous 
individual was recorded as pbest. The best location of the 
herd is the gbest. The herd was randomly generated from 
the population. Finding the best solution by updating the 
velocity and position of each individual according to 
equations (4) and (5). 


x t+1 = x 1 + v t+1 

A id A id ^ v id 


(4) 


4 +1 =w*v‘ id + cfu * (v id - 4 ) + (5) 

+c 2 r 2i *( Pgrf-4) 


Where: t is the t th iteration of the search process, d is the 
size in the search space, deD. ci, C 2 are acceleration 
constants, m, r 2 i are random values, valid in the range 
[0,1]. pid and p g d represent pbest and gbest particles of size 
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d. w is the inertial weight, is the velocity, limited to the 
maximum velocity vmax, vL e[-vmax, vmax]. 


The original PSO algorithm applied to the problem of 
continuity. Kennedy and Eberhart developed the BPSO 
algorithm for the discrete problem (Kennedy & Eberhart, 
1997). The velocity in BPSO represents the element that 
can take the value 1. Equation (4) is still used to update the 
velocity while x lt i, pid get the value 0 or 1. The function 
sigmoid s(v;d) is used to convert the value of Vid into a 
range of values (0,1). The BPSO updates each instance's 
position using equations (6) and (7). 


f 1, if rand() < s(v jd ) 
j 0, otherwise 


( 6 ) 


s(y id ) 


l 

l + £>~ Vi " 


(7) 


The function rand() is a random function whose value is 
distributed in (0,1). 

BPSO 

Begin 

Data set, kfold; D: dimensionality of search space 
N : the population size; T : maximum iterations; 

Cl, C2, Vmax, W 

Randomly initialise the position and velocity of each 
particle; 

while t <T 

evaluate fitness of each particle according to 
Equation (1); 

for i=l to N 

update the pbest of particle i\ 
update the gbest of particle i\ 

end 

for i=l to N 
for d=l to D 

update the velocity of particle i according to 
Equation (4); 

update the position of particle i according to 
Equations (6) and (7); 

end 

end 

end 


return the classification error of the selected feature 
subset; 


end 


IV. RESULTS 
Features and samples: 

The study was tested on the IEEE 39-bus scheme. It 
includes 39 buses, 19 loads, 10 generators. The diagram 
IEEE 39-bus scheme is shown as Fig. 2. It was used in 
many published works. The off-line simulation was 
implemented to collect data for training. Load levels are 
(20,30,...,120)% normal load. The setting fault clearing 
time (FCT) is 50ms (Glover, Sarma, & Overbye, 2012). In 
this paper, all kinds of faults such as single phase to 
ground, double phase to ground, three phases to ground 
and phase-to-phase short-circuit are considered. Faults are 
tested in any buses and in each of 5% distances of long 
transmission lines of the test systems. For each of the 
considered load samples, the generator samples have been 
got accordingly by running optimal power flow (OPF) tool 
of Power-World software. 

The input and output feature are 
x {delVbus.deIPLoad.dcIPn, iw } and y{ 1.0}. Total of input 
features is 104, x{ 104(39+19+46)}. The number of output 
feature is one, y {1.0}. From simulating results, there are 
1167 samples that include 834 S samples and 783 U 
samples. The K-NN (K=l) is used as a classifier for 
evaluating accurate classification because of its simplicity. 
Assessment of classification error is a cross-assessment 
method. It is a 10-fold cross-validation. 



calculate the classification error of the selected feature 
subset; 

return the position of gbest (the selected feature 
subset); 


Fig. 2: The IEEE 39-bus diagram 

Feature selection results: 

The BPSO&l-NN algorithm works with different N 
values, namely 10, 20, 30, 40, and 50. The values of w are 
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0.4, 0.5, 0.6, 0.7, 0.8, and 0.9. The number of iterations is 
100 for the program executions, T=100. The values of ci 
and ci are selected unchanged during program execution, 
Ci=2, C2=2. The program is also executed on Matlab 2018a 
software. Identification tool, 1-NN, is supported from this 
software. 

The convergence characteristics of BPSO&l-NN variable 
selection algorithm are shown as Fig. 3 to Fig. 8. The 
results of the selected number of variables, the value of the 
best fitness function and the test error are presented in 
Table 1. 

Table 1. Best results during the implementation of BPSO 
program with w = 0.9, N = 50, the number of variables in 
the set of variables is 23, the identification error is 5.07%. 
Thus, the algorithm BPSO&1-1NN has reduced the 
number of variables from 104 to 23 variables. The number 
of variables is reduced by 78% while test accuracy reaches 
94.93%. This is the accepted result in previously published 
(Abdelaziz & El-Dessouki, 2013; Amjady & Majedi, 
2007; Dong et al., 2012; Haidar, Mustafa, Ibrahim, & 
Ahmed, 2011; Kalyani & Swarup, 2013; Karami & 
Esmaili, 2013). 


Table 1. Results of BPSO algorithm 


w 

N 

f 

Best Fitness 

ErrsF(%) 

0.4 

10 

32 

7.3868e-04 

5,07 


20 

32 

7.2150e-04 

5,07 


30 

33 

7.3166e-04 

5,01 


40 

30 

6.7693e-04 

5,13 


50 

37 

7.1073e-04 

5,01 

0.5 

10 

33 

7.1424e-04 

5,19 


20 

32 

6.9573e-04 

5,26 


30 

33 

7.0553e-04 

5,07 


40 

34 

7.2445e-04 

5,07 


50 

31 

7.3703e-04 

5,26 

0.6 

10 

33 

7.2295e-04 

4,95 


20 

32 

7.2150e-04 

5,19 


30 

36 

7.3666e-04 

5,01 


40 

33 

6.5327e-04 

4,89 


50 

36 

6.9119e-04 

5,13 

0.7 

10 

31 

7.0315e-04 

5,13 


20 

35 

7.1702e-04 

5,26 


30 

30 

6.9364e-04 

5,26 


40 

29 

7.0089e-04 

5,13 



50 

26 

6.6600e-04 

4,95 

0.8 

10 

32 

6.6138e-04 

5,07 


20 

30 

6.6021e-04 

4,89 


30 

33 

6.5327e-04 

5,26 


40 

30 

6.6857e-04 

5,07 


50 

27 

6.1843e-04 

5,26 

0.9 

10 

32 

6.6996e-04 

5,44 


20 

27 

6.2646e-04 

5,38 


30 

27 

6.1040e-04 

5,13 


40 

23 

5.9552e-04 

5,13 


50 

23 

6.1079e-04 

5,07 




6.5 - 1 - 1 - 1 - 1 - 1 - 1 -'- 1 -'- 

0 20 40 60 80 100 

T=100, w=0.4 


Fig. 3: Convergence characteristics of BPSO&l-NN, 
w=0,4 



Fig. 4: Convergence characteristics of BPSO&l-NN, 
w=0,5 



T=100, w=0.6 

Fig. 5: Convergence characteristics of BPSO&l-NN, 

w=0,6 


www.iiaems.com 


Page | 435 


































































































International Journal of Advanced Engineering, Management and Science (IJAEMS) 
h ttos://dx. doi. ora/10.22161/iiaems. 69 .3 


[Vol-6, Issue-9, Sep-2020] 
ISSN: 2454-1311 



Fig. 6: Convergence characteristics of BPSO&l-NN, 
w=0,7 



Fig. 7: Convergence characteristics of BPSO&l-NN, 
w=0,8 



Fig. 8: Convergence characteristics of BPSO&l-NN, 
w=0,9 

V. CONCLUSION 

The paper successfully applied BPSO algorithm combined 
with 1-NN identifier in the variable selection approach to 
the problem of identifying stability of power system. Test 
results on IEEE-39 bus diagram showed that the number of 
variables decreased significantly, but the test accuracy still 
met expectations. This achieved result will be the driving 
force for the research and application of evolutionary 
algorithms to select variables for the problem of 
identifying stability of power systems. 
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